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ABSTRACT 

We report the first mass and distance measurement of a caustic-crossing bi¬ 
nary system OGLE-2014-BLG-1050L using the space-based microlens parallax 
method. Spitzer captured the second caustic-crossing of the event, which oc¬ 
curred ~10 days before that seen from Earth. Due to the coincidence that the 
source-lens relative motion was almost parallel to the direction of the binary-lens 
axis, the four-fold degeneracy, which was known before only to occur in single- 
lens events, persists in this case, leading to either a lower-mass (0.2 M 0 and 
0.07 Mq) binary at ~1.1 kpc or a higher-mass (0.9 M 0 and 0.35 M 0 ) binary at 
~3.5 kpc. However, the latter solution is strongly preferred for reasons includ¬ 
ing blending and lensing probability. OGLE-2014-BLG-1050L demonstrates the 
power of microlens parallax in probing stellar and substellar binaries. 

Subject headings: gravitational lensing: micro — stars: binary 


1. Introduction 

The detection of binary star systems depends on a combination of diverse observational 
techniques. For example, nearby wide binaries can be directly resolved by high-resolution 
imaging, while close binaries can be detected via eclipsing or spectroscopic methods. To get a 
full picture of the distributions of the mass ratios, masses and separations of binary systems, 
one also needs such a technique as microlensing to probe those binary systems that are 
difficult for other techniques, such as very low-mass binaries (i.e., brown dwarf binaries), dark 
binaries (e.g., binary black holes), and normal binaries with intermediate separations. For 
example, two brown dwarf binaries, OGLE-2009-BLG-151/MOA-2009-BLG-232 and OGLE- 
2011-BLG-0420, were d etected via microlensing, with reported total masses of 0.025 M 0 and 
0.034 M 0 , respectively f Choi et ah 2013 1. 


The challenge faced by standard microlensing observations is to break the degeneracy 
between the mass of and the distance to the lens system, since these two physical parameters 
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both enter a single observable quantity - the timescale of the microlensing event, t^. 

mas 


£e — -; #E — \f K = 


hge 


c 2 AU 


8.14- 


M, 
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Here 6e is the angular Einstein ring radius, /i geo is the geocentric lens-source relative proper 
motion, and 

TTrel = 7T L “ 7T S = AU ( - -=^- 

V^L MS, 

is the lens-source relative parallax. 


The so-called “microlens parallax” can in principle be used to meet this challenge (jGouJd 


19921 ). since the measurement of the microlens-parallax amplitude, 7 Te = 7r re i/#E, directly 
leads to the determination of lens mass Ml and distance D l = AU/Vl, by 


M l = 


e F 


K7Te 


7rL — TrE^E + • 


( 1 ) 


To produce precise constraints on Ml and D^, one therefore needs precise measurements of 
both 6e and tte (keeping in mind that tts is usually known quite accurately). 

Two broad classes of methods have been proposed to measure the microlens parallax 
7 Te. The first is to obtain ob servat ions from a single platform that is being accelerated, 


d be 


which cou 


chron ous (Gould 


(e.g., 


arth (IGouldl 119921) . or a satellite in low-Earth (IHonmal Il999h or geosyn- 
2013h orbit. This method ha s already produced ~100 7 Te measurements 


Alcock et al. 


1995c IPoindexter et al.l 12005c iGaudi et al.l 120081) . However, it is strongly 


b iase d, in the case of binary-lens events, toward nearby lenses flChoi et alJl2013 ; J un g et al. 


20151 ). The ot her class of methods is to obtain obse r vations from at least two well-separated 


observatories (jRefsdallll966l : iHardv & Walker! Il995l : lGouldlll997l ). In order to produce sub¬ 


stantially different light curves observed by different observatories, which lead to precise 
measurement of 7 Te, the required separation between observatories should be ~ AU. For this 
reason, the so-called “terrestrial microlens parallax” 


tori es at different sites, only works in very rare cases ( Gould et al.l 2009 


-ju_—-r? 

(Gould et al. 

2009; 

Yee et al. 

2009 


Gould fe Yeel 2012). The combination of ground-based observations and space-based obser¬ 
vations from a satellite in solar orbit, called “space-based microlens parallax”, has therefore 
been considered the only way to routinely measure the microlens parallax 7 Te for a substantial 
fraction of all microlensing events ( Refsdalll 19661 ). 


There are also several methods proposed to measure the Einstein ring radius 6e■ In 
principle, 6e ca n be measured by resolving the light centroid of the multiple microlensed 
images (IWalkerl 119951) . However, such a method is presently inaccessible given the fact that 
the typical Galactic microlensing event has 6e ~ mas, although a future space telescope 
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with ultra-precise astrometry may be able to achieve that (GouLd & Yee 2 0141 ). The second 
method is to measure the geocentric lens-source relative proper motion, p geo , when the 
lens and source are well separated, typically a few years before or after the microlensin g 
event occurs, since combining this measurement with tg also yields 6 *e (lAlcock et all 2001)). 
However, since this method relies on light from both the lens and source, it does not work for a 
specific class of interesting events, i.e., events with very low-mass objects or stellar remnants 
as lenses. Nevertheless, this app roach has been used several times to measure physica l 


parameters of microlensing events ( Alcock et al. 2001 ; Bennett et ah 2006[ 


and is likely to be used much more frequently in the future fIGouldl 2 0141 ). At present, the 


Dong e_t ah 20 09). 


most widely exploited method of measuring 6 *e is using the finite source effect that is enabled 
when the source crosses or closely passes by the caustic structures of the lens system. 

What makes binary star systems favorable targets for microlens parallax observations 
is that the finite source effect is often detected in binary lens events due to their relatively 
large caustics, producing precise measurement of $e and, if the microlens p arallax is also 
measured, the lens mass Ml and distance (also see i Graff fe Gouldl 120021 ). It is for this 
reason that special attention was paid to binary lens events when we were granted Director’s 
Discretionary Time for a 100 hr pilot program to determine the feasibility of using Spitzer 
as a paral lax sate llite, altho ugh the main objective was to measure lens masses in planetary 
events (Udalski et ah 2 014 ). 


Here we report on the binary lens QGLE-2 014- BLG- 1050L , which is the second space- 
based parallax measurement of a binary lens ( Don g et ah 120071 ) but the first such measure¬ 
ment for a caustic-crossing binary-lens event, or indeed, any caustic-crossing event. We give a 
summary of the observations from the ground and Spitzer in Section [2j The lightcurve mod¬ 
eling is demonstrated in Section [31 the source and blend characterizations and the derivations 
of physical parameters are presented in Section Q] and Section O respectively. In Section EO 
we present a discussion of our results. 


2. Observations 

The source star for microlensing event OGLE-2014-BLG-1050 lies toward the Galactic 
Bulge held with equatorial and Galactic coordinates (RA, Dec) 2 ooo = (17 h 45 m 07!83, —22°54 / 20 , /0) 
and (/, 6)2000 = (5?09,3?23), respectively. It therefore lies just 0?49 above the ecliptic plane. 
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2.1. Ground-based Observations 


The Optical Gravitational Lensing Experiment (OGLE) collaboration alerted the com¬ 
munity to the new microlensing event OGLE-2014-BLG-1050 on 2014 June 6 (HJD' = 
HJD—2450000 = 6815.3), just 35 hours before OGLE detected its first point on the caustic 
entrance at HJD' = 6816.76, based on observations with the 1.4 deg 2 came ra on its 1.3m 
Warsaw Telescope at the Las Campanas Observatory in Chile flUdalski] 20031 ). 


OGLE-2014-BLG-1050 did not attract much attention when it was initially discovered. 
This situation changed after early modeling of online OGLE data revealed it was a binary- 
lens event. Hence, in order to support the Spitzer observations, the Microlensing Follow-Up 
Network (pFUN) started observing this event beginning HJD' = 6822 using the 1.3 meter 
SMARTS telescope at Cerro Tololo Inter American Observatory (CTIO) (see Section l2J2]h 
Observations were taken only about once per night in the following few days. For this 
event, the most crucial task of follow-up teams is to capture the caustic exit with intensive 
observations so as to constrain d E (see Section [T]). However, it is very difficult to predict, 
based on modeling of the collected data, the exact time of caustic crossing, which normally 
lasts only a few hours. To achieve that goal, /iFUN observed at one hour cadence and 
frequently reviewed these data during the night once the lightcurve began its slow rise toward 
the caustic exit. The cadence for OGLE-2014-BLG-1050 was increased to 10 per hour once 
the caustic exit was recognized. 


2.2. Spitzer Observations 


The gene ral description of the Spitzer observations of this pilot program is given in 


Ud alsk i ejt ah 1 2014 1. In short, the program was organized into 2.6 hour windows roughly 


once per day from June 5 (HJD' = 6814) to July 12 (HJD' = 6851); observing targets had to 
be chosen and submitted by J.C.Y. and A.G. to the Spitzer Science Center ~ 3 days before 
the next observing run once per week. 


For the particular case of OGLE-2014-BLG-1050, OGLE had not yet issued its alert 
when the decisions were made for the first week of observations on June 2 (HJD' = 6811.1). 
On HJD' = 6817.9, M.D. suggested OGLE-2014-BLG-1050 might be a binary-lens event, 
but this was too close to the deadline (6818.1) to modify the observing protocol prior to 
upload. During the following week two modelers (M.D. and V.B.) confirmed it as a binary- 
lens event, and therefore high priority was assigned to this event, leading to twice-per-day 
and once-per-day observations in weeks 3 and 4, respectively, and once or twice per day in 
week 5 until it moved beyond Spitzer 's Sun-angle window at HJD' = 6846.56. 
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In total, we obtained 31 observations from Spitzer from HJD' = 6828 to HJD' = 6846, 
which also happened to capture the caustic exit of the event. As shown in Figure [1] the 
Spitzer light curve shows very similar shape as the ground-based light curve, with the only 
significant difference being an offset of ~10 days in time due to the microlens parallax effect, 
since Spitzer was displaced from the Earth by ~1 AU. 


3. Light Curve Modeling 


The observed ground-based and space-based light curves of OGLE-2014-BLG-1050 are 
shown in Figure [U Only one major feature is noticed in the OGLE data: the prominent 
U-shaped trough, indicating a typical caustic-crossing binary-lens event. The Spitzer data, 
on the other hand, only captured the caustic exit. The fact that the two lightcurves show 
similar shape and almost equal amplitude suggests that the trajectory of the source-lens 
relative motion as seen from Spitzer should have similar impact parameter u 0 as that seen 
from Earth. 


The four-fold degeneracy, coming from the fact that the satellite and Earth can pass on 
the same side, which we denote as (+, +) and ( —, —) soluti ons wit h the two signs indicating 
the signs of u 0 as seen from Earth and Spitzer, respectively ( Gouldll2004j h or opposite sides of 
the lens, which we denote as (+, —) and solutions accordingly, was well invest i gated 

in the case of s ingle-lens microlensing (IRefsdall Il966l: iGouldl 11994 I Gould fe Horne 


Yee et al. 2014 : Calchi Novati et al. 2014h . E 


2013 


In the case of b inary l enses in which the 
binary features are detected by both observatories, for example Udalski et al. (2 01411 . this 
four-fold degeneracy normally collapses to the more general two-fold degeneracy, namely u 0 ,± 
degeneracy, w ho se two solutions give very similar amplitudes of 7r E and therefore D L and M L 
flGould fe Horne 20131). However, due to the fact that Spitzer covered only a small portion 
of the caustic-affected lightcurve, the persistence of such a four-fold degeneracy is suspected 
and confirmed in the current binary event. 


Considering the characteristics of this event, we adopt a small variant of the standard 
modeling parameterization. The modeling of the light curves caused by binary lenses usually 
requires 11 system parameters, including seven basic parameters - the time of the closest 
lens-source approach, £ 0 ; the impact parameter normalized by the Einstein radius, u 0 ; the 
Einstein timescale, fg; the normalized source size, p; the normalized projected separation 


1 W e note that these four solutions are defined in a different way in [Yee et al.l (12014 ; Calchi Novati et al 


(1201411 . The four solutions (+, +), (+, —), (—, —) and (—, +) by our definition correspond to (—, +), (+, +), 
(—, —) and (+, —) solutions by their definition, respectively. 

























































7 


between the binary components, s; the mass ratio, q\ and the angle from the binary-lens 
axis to the lens-source relative motion, a - and four parameters to allow for higher-order 
effects, namely the microlens parallax, quantified by tte,n and 7 Tf r , and the bina ry-lens 
orbital motion, quantified by da/dt and ds/dt. We refer the reader to lUdalsk i et al. (120141 ) 
for the sign definitions of parameters uo, a and da/dt. In addition, each observatory requires 
two flux parameters (/s, /b), so that the total flux is determined by 

/tot(f) = /s ■ A(t) + / B , 

where A(t) is the magnification of the source as a function of time t. These flux parameters 
are defined in a system in which f = 1 corresponds to an 18 mag star. For the current event, 
however, noticing that the time of caustic exit (as observed from Earth), t ce , is much better 
determined than t 0 , we replace t 0 with t ce as a free parameter, and t hen sea rch for the t 0 
that produces the given t ce in order to compute the light curve (ICassanl 120081 ). Parameters 
leading to an unbound binary system, quantitatively /3 > 1, are not taken into account, in 
which ( Dong et ah 2009 1 


/3 = 


E] 


kin 


E„ 


KM Q (yi) 


7t e s 3 7 2 


v pot / x 

and the angular Einstein radius 


_ / ds/dt da 

r 2 # E (7r E + 7r s /6» E ) 3 ’ 7 7 II ,7_L V s ’ dt 


( 2 ) 


0k — 


^*,fid(/s//s,fid ) 1/2 


p 


where #*,fid is the fiducial angular source size that is determined for a fiducial source flux 
/s,fid (see Section H] for the final determination of 0*). 

The linear limb darkening coefficients we ado pt, based o n the source color, are (Tj, T 3 . 6 ) = 
(0.43,0.16), following the same procedures as in lUdalski et ah (2 01411 . 


The Markov Chain Monte Carlo (MCMC) method is implemented to fold the minimum 
and the likelih ood distribution of parameters. As usu al, the point source, quadrupole and 
hexadecapole ( Peicha fe Hevrovskv 20091: Gould 2008h approximations are used when the 
source is approaching, although still reasonably far (a few source radii) from, the caustics. 
For epochs that are near or on crossing caustics, we use c ontour integration, in whic h the limb 
darke ning effect is accommodated by using 10 annuli (IGould fe Gaucherell 119971: IDominik 
1998h . @ However, this contour integration may fail at some particular p oints, in which case 
we use the more time-consuming inverse ray shooting (iDong et al.1120061 ). 


2 A more advanced and optimized version of the contour integration proposed by Bozza ( 201C)h has been 


used in the modeling of real time data by V.B.. 

















































We first fit only the ground-based data, i.e., the I band data from OGLE and CTIO, in 
order to obtain a basic sense of this event. The best-fit parameters of the two solutions (w 0 ,±) 
are listed in Tabled] Note in particular that the parallax parameters ke,n and 7 t E jE , are not 
significantly detected. Because the so urce lies very close to the ecliptic, the two solutions 
suffer from the “ecliptic degeneracy” flJiang et al.l 20 041; ISkowron et al. 2 0111) . The Spitzer 
data were then included, which yields much better constraints on the microlens parallax 
vector 7t e = (he,n, ^e,e)- To do so, one needs a careful initial setup of 7 t EjA t and vr E)E . That 
is, the Spitzer caustic crossing clearly occurs between the data points on HJD' = 6839.03 
and HJD' = 6839.94, so if the trial solutions do not have this property the y 2 minimization 
procedure will never arrive there. We have 


AU (At 0 A 


(3) 


in which A t 0 = £ 0iSat — bun- A% = u n . at , — and D± is the projected separation vector of 


the Earth and satellite flUdalski et al.ll2014 ). For solutions in which Earth and Spitzer pass 
the lens on the same side, the similarity between the ground-based and Spitzer lightcurves 
suggests that A u 0 ~ 0, which leads to tte,n ~ 0. Then, since D± ps 1 AU, A to ~ 10 days, 
and i E ~ 80 days, one obtains tte,e ~ 0.13. For the other two solutions, the situation is 
less straightforward because of the difficulty in determining the direction of Dj_ by simple 
inspection. In principle, one can still estimate 7 t e by working out the geometry more care¬ 
fully, but a much easier approach is to conduct a grid search on (vt Ei at, vr EjE ), with all other 
parameters initially set at the best-fit to ground-based data, so as to find a reasonably good 
starting points for MCMC sampling. The best-fit parameters for all four solutions of the fit 
to the combined data sets are listed in Table [2] and the caustic structures and the lens-source 
relative motion as seen from Earth and Spitzer are shown in Figure [2] 

As expected, the inclusion of Spitzer data reduces the error bars on it e ,n and 7 t E jE 
significantly. The four-fold degeneracy does persist in this binary-lens event, and the four 
solutions have nearly equal y 2 , mostly due to the coincidence that the source-lens relative 
motion is almost parallel to the binary-lens direction, i.e., a ~ 180°, and also the fact that 
Spitzer data do not have a long enough time baseline. However, this four-fold degeneracy is 
effectively broken by other considerations, as we will discuss in Section [5] 

The inclusion of the lens orbital motion effect is important in order to uncover the true 
physical parameters of the lens system ( Park et ahl 20131). For OGLE-2014-BLG-1050, the 
orbital motion effect is only marginally detected (Ay 2 ps 5), but it enlarges the uncertainty 
on parameters such as q and u 0 by a factor of two to three. 

It is interesting to compare this binary-lens event with the on ly planetary event so far 
found with the same method, OGLE-2014-BLG-0124 f Udalski et ah 2014 1. One noticeable 
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difference between these two events is the uncertainty in 7 Te- In the present binary-lens event, 
the uncertainty in 7 t e (10 — 20%) after the inclusion of Spitzer data is considerably larger 
than the one for the planetary event OGLE-2014-BLG-0124 (~ 2.5%). The reason is that 
the Spitzer light curve of OGLE-2014-BLG-1050 has fewer features, which itself derives from 
two facts. The first is that due to the Sun-angle limitation, Spitzer 's total time baseline is not 
long enough to cover the caustic entrance as was the case for the planetary event. The second 
is that the current program limits the observation cadence to no more than once per day, and 
thus we were not able to capture the details of the caustic exit. Another consequence of this 
single-feature Spitzer light curve is the asymmetric posterior distributions of 7 Te,at and tie,e 
in the (+,+) and (—, —) solutions. See Figure [3] for the 2D posteriors of (7 t e ,v, tte,®) from 
the (+, +) solution with respect to that from the (+, —) solution. This is because Spitzer 's 
view of the source-lens relative trajectory in the (+,+) [( —, —)] solution can go further in 
the direction of South (North) than North (South). However, this restriction is much weaker 
for the other two solutions (see Figure [2j). We show in the Appendix the triangle diagrams 
of the fitting parameters and derived physical parameters for the two solutions (+, ±). 


4. Color Magnitude Diagram 


We use a varian t of the standard procedure to determine the angular size of the source 
star flYoo et a l. 2004). The standard procedure usually requires magnified images of the event 
taken in both V and / bands. In the case of OGLE-2014-BLG-1050, /iFUN did not take 
any E-band images when the source was substantially magnified. Inste ad, we det ermine the 

V — I color from the I — H color following the procedures introduced bv lYee et al.l (120131) . by 
taking advantage of the fact that the SMARTS camera takes H- band images simultaneously 
with the /-band images. The instrumental I — H color of the source is first determined by 
linear regression of H on the / flux at various magnifications during the event. The V — I 
color is then found by using the color-color relation derived from nearby field stars. Since 
nearby field stars are mostly redder than the source, we have to extrapolate the color-color 
relation blueward to the source position, which introduces a 0.05 mag error in the source 

V — / color. The instrumental /-band baseline flux of the source is determined from the 
modeling (/s,ctio = 0.1612 ± 0.0016). Therefore, the instrumental color and magnitude of 
the source are determined to be {V — /, /)s,ctio = (—0.73,19.98). We then compare the 
source with the centroid of the red clump and find an offset of A(E — /, /)s = (—0.31 , 3.77). 
With the intrinsic centroi d of the red clump (IBensbv et al.l 12013c iNataf et al.l 120131 ) and a 
distance modulus of 14.44 ( Nataf et ah 20131 ). we determine the dereddened source color and 
magnitude to be (V — /,/)s,o = (0.75,18.09), and Mj = 3.65, making the source a turn-off 
star. 
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To determine the source ar 
empirical color-color relations of 

Lgular size, we tl 

Bessell & Brett 

len convert from V — I to V — K using the 
(19881) l 3 l , applv the color/surface-brightness 

relation of 

Kervella et al. 

(2004 

, and finally find 


6+ = 0.80 ± 0.09 pas . 


The error comes from the extrapolation of the color-color relation to the source regime 
0.05 mag), the determination of the red clump centroid using nearby fi eld stars (~ 0.05 
mag) and the derivation of the intrinsic source color (0.05 mag. iBensby et al. 20131 ). 


The lightcurve modeling shows that the source is severely blended. Hence, we also 
determine the color and magnitude of the “blend” (blending light) so as to test whether it 
can be explained by the lens system. From OGLE IV field star photometry, the “baseline 
object” (source plus blend) is found to be (V — /, /)base,OGLE = (2.09,17.81), and the centroid 
of the red clump (V — 7,/)rc,ogle = (2.75,16.27). The correction of non-standard V band 
in OGLE IV is then applied 


A(V — I) j C — 0.92 ■ A(V — -Oogle, 


in which ‘JC’ represents the standard Johnson-Cousins system. With the dereddened red 
clump at (V — 1,1) rc,jc = (1.06,14.32) and assuming the sa me ext inction law Rj = 
AijEiy — /) = 1.23 (the distance modulus and Ry are taken from Nataf et al. 120131 . at the 
event position), the total flux baseline and the source are determined to be (V — /, /)base,jc = 
(2.04,17.81) and ( V — I, /)s,jc = (2.44, 20.08) respectively. We then find the color and mag¬ 
nitude of the blend to be (V — /, /)b,jc = (2-01,17.95). 


We show in Figure [4] the CMD that is used to characterize the source and blend, after 
correction to the standard Johnson-Cousins system, together with the color and magnitude 
of an example lens system from the preferred physical solution (see Section [5]). 


5. Physical Parameters 

The physical parameters derived from our modeling are given in Table HJ in which Uhei = 
(vhei,JV) ^hei,s) is the projected velocity of the lens relative to the source in the heliocentric 


3 In principle, one can directly convert the I — H color to the V — K color. This is not adopted in the 
present work since the centroid of the red clump is not determined in the (/, I — H) plane so well as in the 
(/, V — I) plane, and the determination of the red clump centroid in the (/, / — H ) plane is not within the 
scope of the current work. 
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frame, coming from 


tt e AU 

^hel — -27- V e,J- 


with V 0 j _l ~ (0.7, 28.3) km s _1 being the velocity of Earth projected on the sky at the peak 
of the event. The four-fold degeneracy basically collapses to two physical solutions, since the 
two solutions (+, +) and ( —, —) [collectively (±, ±)] have the same amplitude of parallax 7 t e 
[ as do the other two solutions (±, q=)] and therefore lead to similar lens system properties. 


Our results show that the microlensing event OGLE-2014-BLG-1050 was produced by 
a binary system consisting of 0.9 M e and 0.35 M & stars separated by 5 AU at a distance 
3.5 kpc from the Sun (high-mass binary solution), or consisting of 0.2M 0 and O.O7M 0 stars 
separated by 1.6 AU at a distance 1.1 kpc from the Sun (low-mass binary solution). 


The Rich argument, originally su g gested b y James Rich (circa 1997, private commu¬ 
nication) and elaborated in Calch i Novati et al. (20141). argues that if the two components 
of 7 t e ZUl/AU (namely At 0 /^E and Auo) are small and of the same order, then the (±, ±) 
solutions are strongly favored over the (±, 7 =) solutions because the latter require fine tuning. 
However, because this argument rests on the axial symmetry of the point-lens geometry, it 
cannot be applied to binary lenses in a straightforward manner. 

As seen in Table [31 the two different physical solutions have nearly equal y 2 . However, 
the color and magnitude of the blend determined in Section [I] strongly support the high- 
mass binary solution. We demonstrate this point qualitatively by taking a Sun-like primary 
with a O.4M 0 secondary at 3.2 kpc, all of which are well within the 1-cr error bar o f the 
high-mass binary solution. The primary would have a color V — I = 0.70 ([Ramir ez et al. 


20121 ) and magnitude Mj = 4.15. For the secondary, we take V — I = 2.5 and Mj = 7.6. 
Together with the assumed distance 3.2 kpc, these give the combined dereddened color and 
magnitude (V — /, I) l,o = (0.74,16.64). With the same Ry as we used in Section 0] assumed, 
in order to match the color of the blend, the lens system should suffer an extinction with 
E(V — /) = 1.27 and therefore Aj = Rj ■ E(V — /) = 1.56. This leads to an apparent / band 
magnitude of 18.20 for the lens system. Although it is still 0.25 mag fainter than the total 
blend, we emphasize that this is just a qualitative demonstration of the consistency between 
the lightcurve modeling and photometry. In fact, many factors could be used to explain this 
discrepancy, such as metallicity, stellar evolution stage, or a more massive and closer binary 
system @ . We show in Figure [I] t his example track of the lens system on the CMD. 


4 Note that a more massive and closer binary system is possible since the error bars on the derived mass 
and distance are relatively large. For example, raising (9 E by 5% increases the mass and relative parallax 
both by 5%, leading to a roughly 23% brightness increment. 
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In principle, the blend can be a nearby field star other than the lens. However, this 
scenario is extremely unlikely for this event. Using the OGLE data taken both before and 
during the event, we are able to independently determine the light centroids of the source 
and the “baseline object”. These give an offset of 42 mas, meaning that the source of the 
excess flux must lie within this limit (to be conservative, we use 100 mas for the following 
estimate). Note that this difference (0.15 pixels) is consistent with being coincident because 
there is a star of comparable brightness to the blend that lies just 0.9" away. With the 
measured color and magnitude of the excess flux, (V — I, /)b,jc = (2.01,17.95), we search for 
stars whose brightness is within /b,jc ± 0.5 mag and color within (V — /)b,jc ± 0.3, and find 
344 such stars from the CMD shown in Figure [fl Recalling that this CMD is made using all 
stars within a 3.8' x 3.8' square centered on the microlensing event, this indicates that the 
probability for such a star to lie within 100 mas from the source is only 0.02%. 


A complementary argument that also supports the high-mass binary solution is the 
relative lensing probability. The event rate T = nav as a function of the independent 
physical variables (Ml, -Dl, /x rel ) for single-lens microlensing event is flBatista et al.l 2 0111) 


d 4 T 


Psingie OC dDLdlogMhd 2^ ~ = ^(R,z)(2R E )v rcl f g {M L ) oc D L v(R, z) f (/x rel ) M l ‘ 


(4) 


where v(R, z ) is the number density of stars at position (A, z ) relative to the Galactic center, 
R e = D l 9 e is the Einstein radius at the lens plane, u re i = D L /j, re y is the lens-source relative 
velocity, /(/x rel ) is the two-dimensional probability function for a given source-lens relative 
proper motion /x rel , and S'(Ml) = M^ a is the Galactic stellar mass function in equal bins of 
logM. We take R , z) in the form of 


u(R,z) oc e -R/R* e -\*\/ H , 


and choose R* ~ 3 kpc and H ~ 250 pc; given the Galactic coordinates of this event, 
v(R, z) yields a factor of 1.4 favoring the high-mass binary solution. We choose a = 1 for 
the power-law index of the mass function. For the relative proper motion, /i re i has equal 
amplitude for each solution, b ut its direction results in a factor of ~ 3 in favoring 

the low-mass binary solution flYee et al.l I2014J : ICalchi Novati et al.l I2014T) . For binary-lens 
event, there should be three factors in addition to those in Equation (|4|) . The multiplicity 
frequency, / m , for G-type (primary of the high-mass sol utio n) and M-type (t he primary of 
low-mass solution) dwarfs, is different by a factor of ~ 2 (IDuchene fe Krausil20131 ). The mass 
ratio q results in a factor of unity, since both solutions have the same q. The semi- major axis 


a also leads to a factor of unity, given an Opik law in log a and that a ~ yj3/2a± ( Zhu et al. 


20141 ). With all the above factors as well as the distance factor £)£ considered, this argument 


leads to a conclusion that the high-mass distant binary is more likely to be microlensed than 
the low-mass close binary by a factor of ~ 2. 
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The blend origin strongly prefers the high-mass binary solution, and this solution is 
also preferred by the relative lensing probability. Nevertheless, we suggest that future direct 
imaging of the event with adaptive optic s will help to reach a definitive conclusion about the 
nature of the lens system (IBatista et al.1 20141) . 


6. Discussion 


After 50 years of dreaming, the concept originally proposed by R efsdal] (119661 ) to probe 
the mass function of Galactic astronomical objects without biases toward brightness is fi¬ 
nally underway. With combined observations from Spitzer and ground, we have shown that 
the mass and dista nce of the microlensing planetary events can be well constrained (see 


Udalski et al.l l2014h . and demonstrated the potential of using m icrolens parallax to probe 


the Galactic distribution of planets (see ICalch i Novat i et al. 2 0141) . 


In the present work, we report on a binary-lens event OGLE-2014-BLG-1050 observed 
in our program to demonstrate the power of using space-based microlens parallax to measure 
the mass and distance of binaries. Binary-lens events attract special attention because the 
finite source effect is often detected during caustic crossings and thus leads to well constrained 
Einstein ring radius 0 E . 


Unlike the planetary event QGLE-2014-BL G-0124 , in which the microlens parallax pa¬ 
rameter 7 Te is well determined (2.5% uncertainty, Udalsk i et al. 2014), 7 t e is only constrained 
to within ~ 20% in the present binary-lens event because of its single-feature Spitzer light 
curve. The caustic entrance was not captured by Spitzer for two reasons. First, by the 
time of OGLE’s alert, the entrance had already occurred as seen by Spitzer. Second, given 
Spitzer ’’s Sun-angle limitation it would have been impossible to extend the total time base¬ 
line to capture the caustic entrance even with an earlier alert. However, the uncertainty in 
7 t e could still have been reduced significantly if more Spitzer observations had been obtained 
during the caustic exit. Nevertheless, we emphasize that with Spitzer data the measurement 
of 7r E is quite secure. By contrast, 7r E is not significantly detected if only ground-based data 
are used. 


Another interesting characteristic of OGLE-2014-BLG-1050 is that the four-fold de¬ 
generacy, which usually appears in single-lens events but has not been investigated in the 
binary-lens case, is unexpectedly present. This is mostly due to the coincidence that the 
source-lens relative motion is close to parallel to the binary-lens direction but also to the 
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fact that Spitzer data did not capture the caustic entrance. □ The resulting four degenerate 
solutions are almost equal in y 2 and lead to two very different physical solutions for the 
mass and distance of the lens system: a binary system consisting of 0.9 M & and 0.35 M & 
stars separated by 5 AU at 3.5 kpc (the high-mass binary solution), or a binary consisting 
of 0.2 M q and O.O7M 0 stars separated by 1.6 AU at 1.1 kpc (low-mass binary solution). 

However, this degeneracy is effectively broken when two other factors are considered. 
The color and magnitude of the blend, after the source contribution is subtracted from the 
total baseline, can be well explained by the high-mass binary solution within its 1-cr error 
bars. By contrast, the chance that a random field star is responsible for this blend is only 
0.02%. The lensing probability estimate also favors the high-mass binary solution by a factor 
of two. 


The microlensing event OGLE-2014-BLG-1050 demonstrates the power of microlens 
parallax in measuring ma ss and distance of binaries. Future spac e-based programs wit h 


Spitzer (IGould et ah 


Euclid (Penny ct al. 


2014h and future missions such as Kepler (K2, Gould fe Horne 2013 h 


20131) and the Wild Field InfraRed Survey Telescope flSpergel et ah 


2013 1) can help draw a full picture of the Galactic distribution of binary systems, from brown 


dwarf binaries to binaries involving black holes. 
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5 We remind that this four-fold degeneracy will disappear for programs with long enough time baseline. 
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Table 1. 


Best-fit parameters for groimd-based only fit. 


Parameters 

Uq > 0 (‘+’ solution) 

u 0 < 0 (‘—’ solution) 

y 2 /dof 

457.6/506 

458.4/506 

fee (HJD'—6851) 

0.6318 ± 0.0005 

0.6317 ± 0.0005 

u 0 

0.357 ± 0.014 

0.354 ± 0.017 

f E (days) 

79.2 ± 6.4 

78.0 ± 4.3 

p (io- 4 ) 

6.09 ± 0.64 

6.23 ± 0.50 

7T E,N 

-0.45 ± 0.40 

0.36 ± 0.45 

Ke ,e 

0.31 ± 0.24 

0.40 ± 0.26 

a (deg) 

180.3 ± 2.4 

179.8 ± 2.3 

ds/dt (yr _1 ) 

- 

- 

da/dt (yr _1 ) 

- 

- 

s 

1.108 ± 0.028 

1.115 ± 0.023 
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Table 2. Best-fit parameters for the fit to the combined datasets (no orbital motion). 


Parameters 

(+, +) 

(- -) 

(+, -) 

( - > +) 

X 2 /dof 

489.5/506 

489.6/506 

489.5/506 

490.4/506 

tee (HJD'—6851) 

0.63174 ± 0.00050 

0.63173 ± 0.00050 

0.63157 ± 0.00049 

0.6317 ± 0.00049 

Uq 

0.352 ± 0.015 

0.353 ± 0.015 

0.358 ± 0.014 

0.361 ± 0.016 

t E (days) 

76.7 ± 3.9 

76.9 ± 4.0 

77.5 ± 3.6 

77.9 ± 4.0 

p (io- 4 ) 

6.25 ± 0.47 

6.22 ± 0.47 

6.15 ± 0.45 

6.07 ± 0.47 

tte ,n 

-0.037 ± 0.031 

0.029 ± 0.032 

-0.577 ± 0.056 

0.571 ± 0.057 

tte ,e 

0.115 ± 0.016 

0.117 ± 0.016 

0.089 ± 0.029 

0.098 ± 0.029 

a (deg) 

180.3 ± 1.6 

180.6 ± 1.6 

178.9 ± 1.8 

178.5 ± 1.7 

ds/dt (yr -1 ) 

- 

- 

- 

- 

ds/dt (yr" 1 ) 

- 

- 

- 

- 

s 

1.123 ± 0.020 

1.123 ± 0.020 

1.116 ± 0.019 

1.114 ± 0.021 

Q 

0.382 ± 0.029 

0.380 ± 0.029 

0.383 ± 0.028 

0.383 ± 0.028 


Table 3. Best-fit parameters for the fit to the combined datasets (with orbital motion). 


Parameters 

(+, +) 

(- -) 

(+> ~) 

( _ > +) 

X 2 /dof 

484.7/506 

485.4/506 

483.6/506 

483.5/506 

to (HJD'—6851) 

0.63170 ± 0.00050 

0.63176 ± 0.00048 

0.63166 ± 0.00050 

0.63181 ± 0.00051 

M 0 

0.352 ± 0.036 

0.347 ± 0.030 

0.343 ± 0.034 

0.331 ± 0.033 

t E (days) 

73.3 ± 4.2 

73.3 ± 4.2 

70.0 ± 3.8 

70.3 ± 4.3 

p (io- 4 ) 

6.05 ± 0.64 

6.11 ± 0.59 

6.08 ± 0.65 

6.06 ± 0.58 

ftE,N 

-0.046 ± 0.051 

0.036 ± 0.036 

-0.571 ± 0.063 

0.556 ± 0.066 

tJ'e ,e 

0.111 ± 0.025 

0.115 ± 0.020 

0.081 ± 0.035 

0.098 ± 0.031 

a (deg) 

178.8 ± 2.1 

178.8 ± 2.0 

176.9 ± 2.3 

176.5 ± 2.1 

ds/dt (yr -1 ) 

-0.10 ± 0.40 

-0.19 ± 0.39 

-0.37 ± 0.49 

-0.60 ± 0.52 

da/dt (yr -1 ) 

0.57 ± 0.27 

-0.53 ± 0.27 

0.90 ± 0.53 

-0.77 ± 0.52 

s 

1.113 ± 0.028 

1.113 ± 0.025 

1.105 ± 0.027 

1.097 ± 0.031 

Q 

0.389 ± 0.094 

0.367 ± 0.078 

0.354 ± 0.093 

0.314 ± 0.090 
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Table 4. Physical parameters (ground-based + Spitzer, with orbital motion). 


Parameters (+, +) (-, —) (+, -) (-,+) 


-^primary (Me) 
-^secondary (M e ) 

D h (kpc) 
a± (AU) 

9e (mas) 
v hel:N (km/s) 
hhei ,e (km/s) 


0.91 ± 0.17 
0.355 ± 0.079 
3.46 ± 0.37 
5.02 ± 0.43 
1.34 ± 0.16 
-75 ± 46 
205 ± 47 


0.94 ± 0.17 
0.342 ± 0.092 
3.47 ± 0.19 
5.04 ± 0.37 
1.32 ± 0.14 
58 ± 46 
208 ± 32 


0.21 ± 0.03 
0.073 ± 0.016 
1.12 ± 0.17 
1.64 ± 0.14 
1.33 ± 0.15 
-42 ± 4 
35 ± 2 


0.22 ± 0.04 
0.067 ± 0.013 
1.15 ± 0.17 
1.66 ± 0.14 
1.32 ± 0.13 
42 ± 5 
37 ± 2 
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HJD-2450000 


Fig. 1.— Spitzer’’?, parallax viewpoint of OGLE-2014-BLG-1050, while located ~ 1 AU 
west of Earth. The very broad U-shaped trough in the ground-based light curve indicates 
a typical caustic-crossing binary-lens event. Spitzer saw the same feature but ~ 10 days 
earlier, suggesting that the lens has a projected velocity of ~ 200 km s" 1 due East if Spitzer 
traces the same trajectory as Earth. The inset shows the details of the caustic exit seen from 
Earth at ~ 6851. 
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0.4 0.2 0.0-0.2-0.4 0.4 0.2 0.0-0.2-0.4 

6>/6> e (East) 9 / 6 e (East) 


Fig. 2.— Caustics (grey) and source-lens relative trajectories as seen from Earth (black line) 
and Spitzer (red line); the red and black dots mark the source positions as seen from Earth 
and Spitzer when Spitzer started its observations on HJD' ~ 6828 (the size of the dots here 
do not represent the source size); in each panel, the plus sign marks the position of the 
primary star, and the secondary is ~ 1.1 Einstein units away along the direction of the axis 
of symmetry and on the opposite side of the caustic. In the favored solutions (upper panels), 
Spitzer traces the same trajectory as Earth but ~ 10 days earlier, resulting in a similar light 
curve but displaced by ~ 10 days (see Figure Q]). However, similar lightcurves could also be 
induced if Spitzer passed the lens on the opposite side but with similar impact parameters 
(lower panels). A much larger 7 Te is required in these cases, leading to the noticeable wiggles 
on the source trajectory as seen from Earth. 
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TTe 'TTe ,e 


Fig. 3.— The 2D posteriors between tt^n and tte,e for the (+, +) and (+, —) solutions. 
The contours (representing 1-cr, 2-cr and 3-cr limits) enclose probabilities of 39%, 86% and 
99%, respectively. Compared to the (+, —) solution, the (+, +) solution has much more 
asymmetric contours in the ^e,e) plane. 
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V-I 


Fig. 4.— The color-magnitude diagram of stars in a 3.8' x 3.8' square centered on the source 
star. The positions of the centroid of the red clump (‘RC’), the source star, the blended 
light, and an example lens system taken from the high-mass binary solution before and after 
the extinction are marked. 
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A. Appendix material 
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Fig. 5.— The full triangle diagram of the fitting parameters in the (+, +) solution without 
orbital motion. The contours (representing 1-cr, 2-cr and 3-cr limits) enclose probabilities of 
39%, 86% and 99%, respectively. 
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i ce (HJD' -6851) u o t E (days) P (1CT 4 ) ^e ,n ,e a (deg) s q 


Fig. 6.— Similar to Figure [5] but for the (+, —) solution. 
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iM' 


D h (kpc) 


a x (AU) 


^hel.JV 


w hel,£ 


Fig. 7.— The full triangle diagram of the derived physical parameters in the (+, +) solution 
without orbital motion. Contours have the same meanings as in Figure [51 
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Fig. 8.— Similar to Figure [7] but for the (+, —) solution. 






















































































